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I.  INTRODUCTION 


The  accurate  prediction  or  experimental  determination  of  projectile  aero¬ 
dynamics  is  of  significant  importance  to  the  shell  designer  and  ballistician. 
The  shell  designer  requires  accurate  aerodynamic  data  for  the  overall  develop¬ 
ment  of  new  shell.  The  ballistician  is  concerned  with  the  development  of  aim¬ 
ing  data  and  therefore  relies  heavily  on  accurate  aerodynamic  data.  Experi¬ 
mental  costs  have  skyrocketed  in  recent  years  and  have  contributed  signifi¬ 
cantly  to  overall  system  development  costs.  Computational  techniques  are 
beginning  to  show  promise  as  a  means  to  alleviate  or  at  least  temper  these 
rising  development  costs  by  providing  relatively  low  cost  computer  analysis  of 
new  designs.  As  computer  technology  increases  and  machines  become  faster  with 
larger  memory,  the  use  of  computational  methods  in  design  becomes  more  of  a 
real ity. 

The  means  to  compute  projectile  aerodynamics  for  all  Mach  number  regimes 
covered  by  a  given  projectile  in  its  flight  history  has  been  an  area  of 
research  actively  pursued  by  the  Aerodynamics  Research  Branch  of  the  Ballistic 
Research  Laboratory.  Early  work  had  focused  on  the  supersonic  flight  regime 
and,  in  particular,  on  the  accurate  prediction  of  the  Magnus  force.  The  Magnus 
force,  which  is  very  small  in  magnitude  (on  the  order  of  1/10  the  normal 
force),  is  a  critical  parameter  in  determining  the  dynamic  stability  of  shell. 
The  Magnus  force  is  generated  by  a  spin-induced  distortion  of  the  boundary 
layer;  therefore,  correct  modeling  of  the  viscous/inviscid  interation  is 
critical  for  accurate  computations.  The  work  of  Sturek  et.  al.1  has  shown 
that  accurate  results  in  the  supersonic  regime  can  be  obtained  for  ogive- 
cylinder  projectile  shapes.  This  technique  involved  separate  computations  of 
the  turbulent,  viscous  boundary  layer  and  the  outer  inviscid  flow  field.  As 
the  projectile  shapes  were  generalized  to  include  boattails,  more  sophisti¬ 
cated  computational  techniques  had  to  be  employed.  These  new  methods,  which 
solved  the  thin-layer  Navier  Stokes  equations,  were  successfully  applied  to 
ogi  ve-cyl i nder-boattai  1  shapes  by  Sturek  and  Schiff.2’3  The  solution  of  the 
Navier  -  Stokes  equations  allows  for  the  simultaneous  computation  of  the 
viscous/inviscid  flow  field  and  thus  provided  the  basis  for  good  Magnus  force 
prediction. 


1.  Sturek ,  W.B. ,  et  al.,  "Computations  of  Magnus  Effects  for  a  Yawed, 

Spinning  Body  of  Revolution AIAA  Jourmal .  Vol.  16,  No.  7,  July  1978,  pp. 
687-692. 

2.  Schiff,  L.B.,  and  Sturek,  W.B. ,  "Numerical  Simulation  of  Steady  Supersonic 
Flow  Over  Cone  Ogive-Cylinder-Boattail  Body ,"  AIAA  Paper  No.  80-0066, 
January  14-16,  1980 ;  also,  ARBRL-TR-02363,  U.S.  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  Aberdeen  Proving  Ground,  MD  21005,  September  1981 

(ad  AioeoVo). 

3.  Sturek,  W.B.,  and  Schiff,  L.B.,  "Computations  of  the  Magnus  Effect  for 
Slender  Bodies  in  Supersonic  Flow,  "  AIAA  Jj2MJ21MlL>  Vol.  20,  No.  12, 

December  1982,  pp.  1724-1731 . 
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A  region  of  critical  aerodynamic  behavior  occurs  in  the  transonic  regime, 
0.90  <  M  <  1.15,  where  aerodynamic  coefficients  have  been  found  to  increase  by 
as  much  as  100%.  This  flight  velocity  regime  is  both  experimentally  and  comp¬ 
utationally  difficult.  Thus,  only  a  small  amount  of  experimental  data  are 
available  for  design  studies  and  only  limited  computational  studies  have  been 
made.  An  initial  attempt  to  develop  a  computational  capability  suitable  for 
Magnus  prediction  at  transonic  velocity  was  made  by  Nietubicz,  et.  al4  whereby 
the  thin-layer  Navier  Stokes  computational  technique  was  applied  to  standard 
and  hollow  projectile  shapes  at  zero  degrees  angle  of  attack.  Computational 
results  were  also  obtained  for  non-spinning  projectiles  at  angle  of  attack^ 
however,  limitations  of  computer  resources  (CDC  7600)  became  apparent.  The 
computational  results  obtained  by  Deiwert,6  on  the  ILLIAC  IV,  indicating  the 
complicated  flow  pattern  which  exists  on  a  boattailed  afterbody,  demonstrates 
the  capability  of  the  Navier-Stokes  codes  given  sufficient  computer  resources. 
In  this  paper,  recent  results  are  presented  for  Magnus  force  computations 
using  the  Cray  IS  computer.  Comparisons  are  made  with  earlier  CDC  computer 
results  and  are  further  compared  to  some  limited  experimental  data.  A  brief 
discussion  of  the  numerical  technique  is  included. 


II.  GOVERNING  EQUATIONS  AND  COMPUTATIONAL  TECHNIQUE 

The  Navier-Stokes  equations  solved  here  make  use  of  the  thin-layer 
approximation.7  That  is,  the  viscous  terms  are  neglected  in  both  the  longitu¬ 
dinal  and  circumferential  directions.  The  viscous  terms  are  retained, 
however,  in  a  direction  nearly  normal  to  the  surface  where  large  flow  field 
gradients  exist.  This  formulation  retains  the  momentum  equations  in  all  three 
coordinate  directions.  The  retention  of  the  three  momentum  equations  allows 
for  the  computation  of  separated  flow  and  thus  differs  significantly  from 
boundary  layer  assumptions. 

The  equations  solved  here  are  written  in  a  generalized  coordinate  system. 
This  allows  a  wide  variety  of  body  shapes  to  be  computed  using  the  same  basic 
numerical  technique.  The  notation  for  the  physical  and  transformed  coordinate 


4.  Nietubicz ,  C.J. ,  Pulliam,  T.H.,  and  Steger,  J.L.,  "numerical  Solution  of 
the  Azimuthal-Invariant  Thin-Layer  Navier-StokeA  Equations,  "  AIAA  Jornmal, 
Vol.  18,  No.  12,  December  1980,  pp.  1411-1412 . 

5.  Nietubicz,  C.J. ,  "Navier-Stokes  Computations  for  Conventional  and  Hollow 
Projectile  Shapes  at  Transonic  Velocities,"  AIAA  Paper  No.  81-1262,  AIAA 
14th  Fluid  and  Plasma  Dynamics  Conference,  Palto  Alto,  CA,  1981j  also, 
ARBRL-MR-03184,  U.S.  Army  Ballistic  Research  Laboratory,  ARRADCOM, 

Aberdeen  Proving  Ground,  MD  21005,  July  1982  (AD  A1168s6). 

6.  Deiwert,  G.S.,  "Numerical  Simulation  of  Three  Dimensional  Boattail  After¬ 
body  Flow  Field, "  AIM  Joirrnalj  Vol.  19,  No.  5,  May  1981,  pp.  582-588. 

7.  Baldwin,  B.S.,  and  Lomax,  H.,  "Thin  Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-257,  January  1978. 
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systems  are  shown  in  Figure  1.  The  three-dimensional,  transformed,  thin-layer 
Navier-Stokes  equations,  written  in  nondimensional ,  strong  conservation  law 
form  are8 


3  q  +  3i  +  3  F  +  3G*  =  Re_1arS 
C  n  5  S 


(1) 


The  general  coordinate  transformations  are  defined  as 


5  =  S(x,y,z,t)  -  longitudinal  coordinate 
n  =  n(x,y,z,t)  -  circumferential  coordinate 
C  =  c(x,y,z,t)  -  near  Normal  coordinate 
x  =  t  -  time 


and 


~  p 

pU 

pu 

-  _i 

E  =  J 

pull  +  5XP 

pv 

PVU  + 

pw 

pwU  +  SZP 

e 

(e+p)U  -  StP 

L  u  - 

pV 

Pw 

puV  +  nxP 

*  -1 

G  =  J  1 

puW  +  ^xp 

pW  +  n  p 

pvW  +  Cyp 

J 

pw V  +  nzP 

' 

pwW  +  ^ZP 

(e+p)V  -  ntP 

(e+p)W  -  s.p 

L  L  -J 

The  viscous  matrix,  S,  is  written  as 


8.  Pulliam,  T.H.,  and  Stegev,  J.L.,  "On  Implicit  Finite-Difference 

Simulations  of  Three-Dimensional  Flow,"  AIAA  Jbirrnal,  Vol .  18,  Do.  2. 
February  1980,  pp.  159-167. 


9 


0 


uU2n2+S2)u  +(u/3)(c  u  +  C  v  H w  )C 


x  'y  'z'  ? 


x  c  *y  c  z  C'  x 


S  =  J 


-1 


{ ( cJ+^+tpCu/2(u2+v2+w2) £+kP r” 1  (y-1  )“ x(a  ) ?] 
+  (u/3)  (^xu+?yv+?zw)  Ux»Vcy  Vczwc^  ^ 


The  velocities 


u  =  «t  +  v +  y +  y 

V  =  nt  +  nxu  +  nyv  +  nzw 

W  =  ct  +  cxu  +  cyv  +  czw 


(2) 


represent  the  contravarient  velocities.  The  nondimensi onal  velocities  U,  V, 
and  W  are  those  components  in  the  direction  of  the  transformed  coordinates 
5,  n,  and  t,  respectively.  The  Cartesian  velocity  components  u,  v,  w  together 
with  the  density  p  and  total  energy  per  unit  volume  e  are  retained  as  the 
dependent  variables.  The  local  pressure,  p,  is  determined  using  the  relation 


p  =  (y  -  1) (e  -  .5p(u2  +  v2  +  w2)). 


The  velocities  are  nondimensionalized  by  the  free-stream  speed  of  sound  a^, 
the  density  by  and  the  total  energy  by  p^a2.  The  additional  parameters 

appearing  in  equation  1  are:  (a)  coefficient  of  thermal  conductivity,  k; 
(b)  dynamic  viscosity,  y;  (c)  Reynolds  number  based  on  body  diameter,  Re; 
(d)  Prandtl  number,  Pr. 

As  mentioned  earlier,  these  equations  are  written  in  transformed  coordi¬ 
nates;  therefore,  the  various  body  shapes  are  introduced  through  determination 
of  the  metric  terms  £  ,  n  ,  C7,  etc.  These  terms  are  formed  by  a  combination 

of  the  derivative  terms  x^,  y^,  z^,  etc.,  and,  together  with  the  transforma¬ 
tion  Jacobian,  allow  for  variable  body  geometries.  Thus,  .one  of  the  first 
steps  in  performing  a  computation  is  the  generation  of  a  computational  grid 
which  provides  the  x,  y,  z  points  for  the  metric  determination.  These  points 
are  determined  prior  to  the  computations  and  are  not  changed  with  time. 
Examples  of  the  computational  grid  used  in  this  study  are  shown  in  Figures  2 
and  3.  A  two-dimensional  slice  of  the  overall  grid  is  shown  in  Figure  2.  The 
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upstream,  downstream  and  outer  flow  field  computational  boundaries  extended 
approximately  18  body  diameters  from  the  body  surface.  At  this  distance  the 
flow  field  should  be  uniform  and  the  imposed  boundary  conditions  are  consi¬ 
dered  valid.  Figure  3  shows  an  expanded  view  near  the  body.  The  clust.  r ing 
of  qrid  points  near  the  body  surface  is  required  in  order  to  resolve  trie 
viscous  components  of  the  flow  field  near  the  body  surface.  Due  to  the  lack 
of  sufficient  computer  storage,  judicious  use  must  be  made  of  the  limited  grid 
points  available.  In  regions  where  the  viscous  effects  are  not  predominant 
and  the  flow  field  changes  slowly,  the  grid  points  are  sparse.  Additional 
grid  clustering  is  used  in  the  longitudinal  direction  where  flow  field  gra- 
dients  are  expected.  The  two-dimensional  grid  shown  in  Figure  2  was  rotated 
about  the  body  axis  in  10 -degree  increments  for  the  Cray  computations  (20- 
deqree  increments  for  the  CDC  7600)  in  order  to  obtain  the  three-dimensional 
grid  required  for  computations  at  angle  of  attack.  The  wake  was  modeled  by 
extending  the  boattail  for  approximately  two  calibers  (X/D  -  7.8),  turning  the 
grid  line  parallel  to  the  model  axis  and  extending  it  to  the  downstream 

boundary. 

As  mentioned  in  the  introduction,  the  Magnus  effect  is  produced  by  a  spin 
induced  distortion  of  the  boundary  layer.  The  computation  must  therefore  e 
fully  three-dimensional  since  no  plane  of  symmetry  exists.  The  boundary 
conditions  used  for  the  computations  are: 

(1)  inner  boundary,  body  surface 

U  =  U  =  0 
V  =  “!L 

■  aoo  ,  . 

p  -  first  order  extrapolation 

p  -  calculated  using  the  three  transformed  momentum 
equations 

(2)  outer  boundary 

constant  free  stream  values  are  used  for  all  variables 


(3)  downstream  boundary 

M  <  1  pressure  is  fixed  at  pM  and  all  other  variables 
are  extrapolated 

M  >  1  first  order  extrapolation  on  all  variables 

The  numerical  scheme  used  for  the  solution  of  Equation  1  is  a  fully 
implicit,  approximately  factored,  finite  difference  algorithm  in  delta  form  as 
analyzed  by  Beam  and  Warming.9  This  scheme  can  be  first  or  second  order 
accurate  in  time  and  second  or  fourth  order  accurate  in  the  three  spatial 


9.  Beam,  R.,  and.  Warming,  R.F. ,  "An  Implicit  Factored  Scheme  for  the 
Compressible  Wavier- Stokes  Equations,"  AIAA  Journal,  Vol.  16,  Wo.  4, 
April  1978,  pp.  393-402. 
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.p  on  the  leeward  ray  is  shown  as 


4.5  x  106.  The  Magnus  measurements  were  compared  f°r samaJJ°^J1m^y0d1ameter 
parameter,  PD/V  =  0.39  where  P  -  mu  comparisons  are  made  using  the 

Langley  dafa”,  while  "Adynamic  coefficient  comparisons  are  made  with  the 
NSWC  data.12 

A.  Surface  Pressure  Coefficient. 

The  surface  pressure  coefficient,  Cp 
function  of  longitudinal I  position  In .  Flgjr.^. 

Tcr CDC  ^r0ensaultrseS(iashe5 Mine)  show  marginal  agreement  with  the  experimental 

velocities. 

Comparisons  between  computation ^1(f«K^t  for^th.  clnamf.rjntl.1 

nM-ainPd5  Dreviously  for  an  axisymmetric  computation,  figure  »  5  yl*  .  n  nf 

Elf  for  thfthree-dimensional  cases  when  adequate  grid  resolution  is  achieved. 

A  primary  purpose  of  this  research  effort  is  the  development t  of fa  capa¬ 
bility  for  the  prediction  of  aerodynamic  coefficients  and,  in  particular, 


12.  Kayser,  L.D..  Ballistic  Research  Laboratory/ARRABCOM,  Aberdeen  Proving 
Ground,  Maryland  21005,  private  aommumeatvone . 
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capability  to  compute  the  Magnus  effect.  As  noted,  the  Magnus  effect  is  a 
viscous  phenonena  associated  with  the  spinning  projectile.  Therefore,  in 
order  for  a  computational  technique  to  predict  this  effect,  it  must  adequately 
compute  the  longitudinal  and  circumferential  wall  shear  stress  for  the  spin¬ 
ning  projectile  at  angle  of  attack.  The  experimental  determination  of  the  u, 
v,  and  w  velocity  distribution  is  especially  difficult  at  transonic  veloci¬ 
ties.  Although  no  experimental  data  are  available  for  comparison,  the  compu¬ 
ted  ci rcumferenti al  velocity  distributions  are  shown  in  Figures  10  and  11  for 
X/D  =  4.22  and  X/D  =  5.50,  respectively.  A  significant  asymmetry  can  be  seen 
in  the  velocity  distributions  at  <j>  =  90°  and  270°.  At  <(>  =  90°  cross-flow 
velocity  caused  by  the  angle  of  attack  is  in  the  same  direction  as  the  wall 
velocity.  At  <j>  =  270°,  the  outer  cross-flow  velocity  opposes  the  wall  velo¬ 
city.  The  ci  rcumferential  velocity  of  profiles  at  4>  =  0°  and  <j>  =  180°  are 
equally  affected  by  the  surface  spin.  Figure  11  shows  the  velocity  distribu¬ 
tion  at  a  station  midway  on  the  boattail.  The  circumferential  velocity  dis¬ 
tribution  at  4>  =  90°  and  <f>  =  270°  has  changed  significantly  from  that  shown  in 
Figure  10.  On  the  boattail,  the  decreasing  body  diameter  results  in  the  sur¬ 
face  velocity  decreasing  in  magnitude.  However,  the  boundary-1 ayer  thickness 
in  this  region  increases  and  the  effect  of  surface  spin  is  seen  to  persist 
further  out. 

B.  Pitch  Plane  and  Magnus  Coefficients 

An  integration  of  the  pressure  and  viscous  forces  has  been  carried  out  in 
order  to  determine  the  aerodynamic  coefficients.  The  sign  convention  used  for 
the  coefficients  is  shown  in  Figure  12.  The  results  in  Figures  13-17  are* 
plotted  as  a  function  of  longitudinal  position  and  thus  show  how  the  force 
develops  over  the  length  of  the  projectile.  Figure  13  is  a  plot  of  the  normal 
force  coefficient  and  shows  the  rapid  increase  in  normal  force  which  occurs  on 
the  ogive  portion  of  the  projectile.  The  cylinder  portion  should  produce  no 
significant  additional  normal  force;  however,  the  computation  indicates  a 
slight  increase  in  normal  force  here.  The  reversal  in  the  direction  of  the 
force  on  the  boattail  can  be  clearly  seen  as  the  accumulated  normal  force 
decreases  over  the  length  of  the  boattail.  The  experimental  normal  force 
coefficient,  indicated  by  the  circle,  shows  very  good  agreement  with  the 
computation.  The  spin  rate  of  typical  artillery  shell  is  of  the  order  of  300- 
500  rev/sec.  As  mentioned  previously,  the  Magnus  effect  results  from  a  spin- 
induced  distortion  of  the  viscous  boundary  laver  which  occurs  for  artillery 
shell  at  angle  of  attack.  Previous  studies*’3  have  shown  that  the  Magnus 
effect  consists  of  the  sum  of  the  boundary -layer  displacement  effect  (asym¬ 
metric  surface  pressure  distribution)  plus  the  viscous  wall  shear  stress  con¬ 
tributions  tx  =  u(-g^-)y_0  and  =  p(^y)y=0-  development  of  all  three  com¬ 

ponents  of  the  Magnus  force  are  shown  in  Figures  14,  15,  and  16,  respectively, 
as  a  function  of  longitudinal  position.  Both  the  longitudinal  and  circumfer¬ 
ential  components  (Figures  14  and  15)  are  seen  to  be  of  the  order  of  10“  5 
and  10"'+,  respectively.  The  pressure  component  (Figure  16)  is  of  the  order 
10”2.  The  dominant  component  of  the  transonic  Magnus  effect  is,  therefore, 

seen  to  be  the  boundary- layer  displacement  effect,  Cyn  .  Additionally,  the 

fw 

largest  portion  of  the  total  Magnus  effect  is  seen  to  develop  on  the  boattail 
where  the  viscous  boundary  layer  reaches  its  maximum  thickness.  This  is  the 
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same  qualitative  behavior  reported  by  Sturek,  et  al1’3  for  supersonic  flow. 

The  total  Magnus  force,  Cy,  and  its  components  Cyu  (longitudinal  shear 

stress),  Cyw  (circumferential  shear  stress),  and  Cyp^  (pressure)  are  shown  in 

Figure  17.  The  total  Magnus  force  can  be  compared  with  the  experimental  data 
whereas  the  components  are  shown  to  indicate  the  relative  magnitude  of  each. 

Considering  the  small  magnitude  of  the  Magnus  force  and  the  agreement  achieved 
for  the  normal  force,  the  quantitative  agreement  between  the  computation  and 
experiment  is  regarded  as  very  good.  The  experimental  Magnus  force  measure-  £ 

ments  were  obtained  in  a  wind  tunnel  not  specifically  designed  for  transonic 
flow  and  are  considered  to  be  of  good  qualitative  value  but  of  questionable 
quantitative  value.  Additional  computations  at  various  transonic  Mach  numbers 
and  reliable  experimental  data  are  required  before  a  full  assessment  of  the 
computational  technique  can  be  made.  This  first  result,  however,  for  predic¬ 
ting  the  Magnus  effect  at  transonic  velocity  is  considered  very  encouraging. 


VII.  SUMMARY 

The  research  effort  presented  in  this  paper  is  part  of  an  overall  program 
to  develop  a  sophisticated  predictive  capability  for  projectile  aerodynamics. 
The  pacing  requirement  for  this  capability  is  the  determination  of  the  Magnus 
force  in  the  transonic  flight  regime. 

An  implicit  finite  difference  code,  which  solves  the  unsteady  thin-layer 
Navier-Stokes  equations,  has  been  applied  to  a  projectile  shape  at  a  =  2.0°,  M 
=  0.91.  The  solution  was  marched  in  time  until  a  steady- state  result  was 
obtained.  Computations  were  first  performed  on  a  CDC  7600  using  a  finite  dif¬ 
ference  grid  of  21,600  points  and  required  7.78  hours  of  computer  time. 
Increased  grid  resolution  with  faster  computational  speed  per  grid  point  was 
obtained  by  performing  the  computation  on  a  Cray  IS  vector  computer. 

The  computations  have  been  compared  to  experimental  surface  pressures  and 
aerodynamic  force  coefficients.  The  circumferential  velocity  distribution, 
presented  for  two  axial  locations,  showed  the  significant  interaction  between 
the  cross-flow  velocity  resulting  from  angle  of  attack  and  the  body-surface 
velocity.  Experimental  velocity  profile  data,  which  are  very  difficult  to 
obtain  for  a  spinning  model  at  transonic  speeds,  are  required  to  fully  assess 
the  computational  results.  The  normal  and  Magnus  fc*rce  coefficients  have  been 
shown  to  be  in  good  quantitative  agreement  with  experimental  data.  The 
individual  components  of  the  Magnus  force  have  additionally  been  presented  and 
indicate  qualitatively  good  results.  The  need  for  additional  grid  resolution 
or  adaptive  grid  techniques13  have  been  identified  as  a  further  requirement  to 
achieve  more  accurate  predictions.  Good  quality  experimental,  transonic 
Magnus  data  is  also  required  for  future  code  validation. 


13.  Dwyer,  H.A. ,  Kee ,  R.J.,  and  Sanders ,  B.R.,  "Adaptive  Grid  Method  for 
Problems  in  Fluid  Mechanics  and  Heat  Transfer A IAA  Journal,  Vol.  18, 
October  1980,  pp.  1205-1212 . 
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The  present  results  indicate  that  the  thin-layer  Navier-Stokes  computa¬ 
tional  technique,  in  conjunction  with  enhanced  computer  technology,  has  the 
potential  of  providing  the  capability  to  accurately  predict  the  aerodynamic 
behavior  of  spinning  shell  at  transonic  velocities,  including  the  Magnus 
effect. 
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Figure  1.  Physical  and  Transformed  Coordinate  System 
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Figure  2.  Physical  Grid  -  Total  Flow  Field 
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Figure  4.  Model  Details 
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Cp,  PRESSURE  COEFFICIENT 


Figure  5.  Surface  Pressure  Coefficient,  Experiment  and 
Computations,  M  =  0.91,  a  =  2.0° 


X/D 

Figure  6.  Windward  and  Leeward  Surface  Pressure 
Coefficient,  M  =  0.91,  a  =  2.0° 
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Figure  7.  Circumferential  Pressure  Distribution, 
M  =  0.91,  a  =  2.0°,  X/D  =  1.56 
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Figure  8.  Circumferential  Pressure  Distribution, 
M  =  0.91,  a  =  2.0°,  X/D  =  5.19 
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Figure  9.  Surface  Pressure  Coefficient,  Experiment  and 
Axisymmetric  Computation,  M  =  0.96,  a  =  0.0° 


t 

Fiqure  10.  Circumferential  Velocity  Profiles  for  <J>  =  0°,  90°,  180°, 
270°  at  M  =  0.91,  a  =  2.0°,  X/D  =  4.22,  PD /V  =  .39 
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Figure  11.  Circumferential  Velocity  Profiles  for  <|>  =  0°,  90°,  180°, 
270°  at  M  =  0.91,  a  =  2.0°,  X/D  =  5.5,  PD/V  =  .39 
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NORMAL  FORCE 


Figure  12.  Aerodynamic  Coefficient  Sign  Convention 
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NORMAL  FORCE  COEFFICIENT 


X/D 


Figure  13.  Normal  Force  Coefficient  Along  the  Projectile, 
Computation  and  Experiment,  M  =  0.91,  a  =  2.0° 


Figure  14.  Circumferential  Wall  Shear  Contribution  to  Magnus 
Force,  M  =  0.91,  a  =  2.0°,  PD /V  =  .39 
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Figure  15.  Longitudinal  Wall  Shear  Contribution  to  the  Magnus 
Force,  M  =  0.91,  a  =  2.0°,  PD/V  =  .39 


Figure  16.  Pressure  Contribution  to  the  Magnus  Force, 
M  =  0.91,  a  =  2.0°,  PD/V  =  .39 
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Figure  17.  Magnus  Components  and  Total  Magnus  Force  Along  Projectile, 
Computation  and  Experiment,  M  =  0.91,  a  =  2.0°,  PD/V  =  .39 
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LIST  OF  SYMBOLS 


a 

speed  of  sound 

A 

reference  area,  ttD 2/4 

CP 

specific  heat  (at  constant  pressure) 

CP 

«* 

CN 

pressure  coefficient 

normal  force  coefficient,  //  Pw  cos<i>  cos0g  dS/qA 

^  Cy 

Magnus  force  coefficient,  Cv  +  Cv  +  Cv 

yu  yw  ypw 

cyp 

//  Pw  si  n<f>  cosOg  dS/qA 

Cy 

•Tu 

Cy 

D 

//  t  si  n<J)  sin6g  dS/qA 
//  cos<t>  cos0g  dS/qA 

body  diameter  (57.15mm) 

e 

total  energy  per  unit  volume/p^a^ 

E,F,G,q 

flux  vector  of  transformed  Navier-Stokes  equations 

J 

transformation  Jacobian 

M 

Mach  number 

P 

pressure/p^a^ 

Pr 

Prandtl  number,  p  c  /« 

7  oo  n  oo 

q 

free-stream  dynamic  pressure 

R 

body  radius 

Re 

Reynolds  number,  pj^a^/p^ 

ReD 

Reynolds  number,  P^Du^/p^ 

•  s 

surface  area 

A 

s 

t 

viscous  flux  vector 

physical  time 

u,v,w 

Cartesian  velocity  components/a^ 

u,v,w 

Contravariant  velocity  components/a^ 
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LIST  OF  SYMBOLS  (cont'd) 


x,y,z  physical  Cartesian  coordinates 
a  angle  of  attack 

Y  ratio  of  specific  heats 

k  coefficient  of  thermal  conductivity/^ 

1 1  coefficient  of  viscosity/y^ 

£,n,c  transformation  coordinates  in  axial,  circumferential  and  radial 

di recti  on s 

p  density/ p,,, 

eD  local  body  slope 

D 

t  transformed  time 

xx  longitudinal  wall  shear  stress 

t  circumferential  wall  shear  stress 

♦ 

<j>  circumferential  angle 

a)  angular  velocity,  rps 

Superscript 

*  critical  value 

Subscri pt 

«  free-stream  conditions 
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